Mean-Field-Theory for Polymers in Mixed Solvents Thermodynamic and Structural 

Properties 



Amina Negadi 1,2 , Anne Sans-Pennincks 3 , M. Benmouna 1,2 , and Thomas A. Vilgis* 1 ' 4 
1 Max-Planck Institut fur Polymerforschung, Ackerrnannweg 10, D-55122 Mainz, Germany 
2 University Aboubakr Belkaid of Tlemcen, Institut of chemistry, BP 119, Tlemcen 13000, Algeria 

3 Centre d'etudes de Bruyeres le Chatel, B.P. 12; F-91680 Bruyeres les Chatel, France 

4 Laboratoire Europeen Associe, ICS, 6, rue Boussingault, F-67083 Strasbourg, France 

February 1, 2008 

Theoretical aspects of polymers in mixed solvents are considered using the Edwards Hamiltonian 
formalism. Thermodynamic and structural properties are investigated and some predictions are 
made when the mixed solvent approaches criticality. Both the single and the many chain problems 
are examined. When the pure mixed solvent is near criticality, addition of a small amount of polymers 
shifts the criticality towards either enhanced compatibility or induced phase separation depending 
upon the value of the parameter describing the interaction asymmetry of the solvents with respect to 
the polymer. The polymer-solvent effective interaction parameter increases strongly when the solvent 
mixture approaches criticality. Accordingly, the apparent excluded volume parameter decreases and 
may vanish or even become negative. Consequently, the polymer undergoes a phase transition 
from a swollen state to an unperturbed state or even take a collapsed configuration. The effective 
potential acting on a test chain in strong solutions is calculated and the concept of Edwards screening 
discussed. Structural properties of ternary mixtures of polymers in mixed solvents are investigated 
within the Edwards Hamiltonian model. It is shown that the effective potential on a test chain in 
strong solutions could be written as an infinite series expansion of terms describing interactions via 
one chain, two chains etc. This summation can be performed following a similar scheme as in the 
Ornstein-Zernike series expansion. 



1. Introduction 

Polymers in mixed solvents, stow different properties than in individual solvents raising interesting questions from a 
fundamental point of viewuO. Specific features of these systems originate from incompatibility of the two solvents 
and asymmetry of their interactions with the polymer. If the system exhibits a preferential adsorption ^Lane solvent 
with the polymer, a difference of solvent composition inside and outside the polymer domain is obs^vedErO. Another 
peculiar behavior of dilute solutions of polymers in binary mixtures of bad solvents is cosolvencyOTla This behavior 
is quite unexpected since two non solvents for the polymer turn out to present together the properties of a good 
solvent which is unusual from the point of view of excluded volume effects. If no specific interactions exist between 
particles in the mixture such as hydrogen bonds or long range electrostatic interaction, one would expect that a 
mixed solvent exhibits a smooth transition from one solvent quality to the other as the solvent composition changes 
from to 1 and the solvent power should be an averaged value of its individual components. This is not the case in 
cosolvent systems where sharp changes characteristic of first order transitions take place. These transitions can be 
monitored through variations of radii of gyration, second virial coefficients, solution viscosities ... with the solvent 
composition. 

In recent years, with the development of new techniques of polymer characterization using molecular labeling, dyes 
and neutron scattering techniques, ...widespread research of polymers in mixed solvents has been witnessed in an 
attempt to understand better their fundamental properties. These mixtures were studied by many authors using 
different methods such as light and neutron scattering. Mixtures of ordinary and deuterated solvents such as light 
and heavy water, ordinary and deuterated organic solvents (benzene, toluene ...) are used quite often in neutron 
scattering measurements to achieve maximum contrast taking advantage of the large difference between the coherent 
scattering lengths of hydrogen and deuterium atomsEJj. 

In polyelectrolyte solutions, a low molecular weight salt is added to increase the ionic strength which weakens the 



electrostatic forces and lowers their range. Salt introduces two additional low molecular weight components (co-ions 
and counter-ions) which could present the properties of mixed solventsLj. r~i r~i n 

Recently, mixed solvents were used to study the thermodynamics of polymer brushesli§c3. Auroy and AuvraytSI 
observed the collapse-stretching transition of grafted polydimcthylsiloxane chains by small angle neutron scattering. 
Changing temperature and composition, they were able to control the height of the brush. The special boundary 
conditions in the terminal attachments result into polymer conformation and dynamics which are different from those 
of the-bulk systems. The degree of swelling of polymer brushes is largely determined by the solvent quality. Rothers 
et al.E3 investigated the swelling of polystyrene (PS) brushes in toluene-methanol mixtures. By changing the solvent 
composition, they were able to see the collapse-stretching transition in the brush. 

Problems involving polymers in mixed solvents are not only of fundamental interest but also they raise interesting 
concrete questions and offer a large spectrum of practical applications. The increase of amount of water into polymer 
domains in nylon technology can beanaoitored directly by NMR imaging techniques and give direct profiles of solvent 
concentration and ingress dynamicsEjO. Research efforts are engaged by polymer industries and rayon manufacture 
to develop new clean methods to meet the requirements of environmental regulations^. In these efforts applications 
of mixed solvent systems are often considered. Polymer synthesis-in mixed solvents is another application which is 
found to allow for a good control of the polymerization processE-3'E£l. 

In this paper, we are interested in the properties of polymers in mixed solvents from thoiundamental point of view 
and present a theoretical investigation based upon the Edwards Hamiltonian formalismc3. In the next section, we 
review the case of a single chain in mixed solvents. The effective potential acting on the chain is calculated together 
with the scattering functions. Phase behavior and chain conformation are obtained from the generalized partition 
function in the Edwards Hamiltonian formalism. In section 3, we consider the case of strong solutions and discuss 
changes in the phase behavior and the scattering properties due to the polymer concentration of the medium. Section 
4 gives concluding remarks. 



2. The single chain in mixed solvents 



A. The Edwards Hamiltonian formalism 



In this section, we consider infinite dilute solutions of polymers in mixed solvents and use the Edwards Hamiltonian 
formalism to determine the effective potential exerted on a chain under various conditions of temperature and solvent 
composition. For simplicity, we assume that all the chains have the same degree of polymerization N and the same 
length I to avoid complications due to polydispersity. For this idealized system, the Edwards Hamiltonian is 



m = ^l ^f) d - + \l I VppKRM -*(»')]<*»<*>' 

pN pN 

+ Y, V AP [n(s)-rf}d S + J2 V BP [R(s)-rf]ds 
i J o i Ja 

1 N 1 

+ ± E ^A[rf - rf] + t J2 VBB [rf rf] + £ V AB [rf rf] (1) 

i,0 i,i hj 

where (3 = 1//cbT, ks is the Boltzmann constant and T the absolute temperature, R(s) is the chain variable and s the 
curvilinear coordinates along the chain ( contour variable). The first term in the right hand side (RHS) of equation 
(1) is the Wiener measure for random walks (assuming that unperturbated chains are Gaussian) and represents the 
chain entropy due to its elasticity. The other terms in the RHS of equation 1 result from two-body interactions 
between various particles in the medium either of the same or different species. In this expression, we have used the 
assumption that the interaction potentials are short ranged and can be approximated with Dirac delta functions 

V i j(T i -T J ) = V ij 5(T i -T j ) (2) 

where Vij is the strength of the interaction potential between particles i and j. Using the standard procedures of 
functional theory and Edwards Hamiltonian formalism, we derive the thermodynamic and the structural properties of 
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the mixture under consideration from the generalized partition function. Starting from the Hamiltonian of equation 
1, one obtains the partition function as a functional integral over all monomers and solvent coordinates 

Z = J VR{s) J JJdr^IJdrf exp -/3W[(R(s), rf , rf ] (3) 

i i 

Where the superscripts A and B refer to the solvents. It is quite difficult to handle the partition function in the 
particle position variable. Rather, it is more convenient to transform into collective variables such as particle densities 

<4 = cxp(-iqif) 

i 

^ = £>p(-iqrP) (4) 

i 

Considering the single chain in solution and integrating over solvent coordinates, one obtains the Edwards Hamilto- 
nian for the chain subject to the effective potential Vb(q) 

m = i?Jo (^) 2 ^+E^ ds / w ^(q) cx p[- if i( R ( s )- R ( s '))] (5) 

The explicit form of Vb(q) is obtained following the procedure described below. 
B. The effective potential of a single chain Vb(q) 

The partition function in equation (3) can be written differently by making use of the collective coordinates and 
Pi 

Z = J VR(s) J n^n^^t-J* jf (^) ds ~\ V ™[ dsds'S[K(s)-K(s')} 

-^apV / cxp(-iqR(s))4-V B pV / exp(-iqR(s))p^ 
„ Jo „ Jo 



q q 



~Vaa £ K\ 2 \v m £ |p«| 2 - V AB Al (6) 



q q q 

where the letter P refers to the polymer. The particle conservation can be expressed by using the following equation 



Pq-^expf-iqr;) 



1 (7) 



q 

and introducing field variables via the auxiliary Edwards random fields <J> 

%q - Yl CX P(~ ic l r i)] = J 5 ®1 CX P R $ qP-q] CX P h i$ q _iqr i] ( 8 ) 

i i 

With straightforward mathematical manipulations, one arrives at 

Z = JvK(s) /lI^qIl^q ex P[-|2 [ (^f) 2 ds-\v PP j" j" dsds>6[R( S )-K(s')] 

r N r N 
-Vap^2 J dscxp(-iqR(s))p£-V B p^ J dsexp(-iqR(s))p^ 



q q 



\ 5> AA + ^)l4| 2 - \ £(Vbb + ^)|^| 2 - Fab E^-ql ( 9 ) 

q 
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For an incompressible mixture, the total mean density fluctuation is zero and one has 



/ dscxp(-iqR(s)) = (10) 
Jo 



Substituting this result into equation (9) and integrating over solvent collective coordinates, one obtains the partition 
function in terms of an effective Hamiltonian TL given by 

0H= wj" i^f) ds + £ £ Vo(q) ex P[" ic l( R ( s ) - WWsds' (11) 
where the effective potential Vb(q) is function of the potentials Vij and the partial structure factors SA(q) and 5b (q) 

i ( Vbb — Vab + Vap — Vbp + q \s I 

<^b (q) Vaa + Vbb - 2 Vab + sKq) + s^) 

The subscripts of potentials run over polymer P and solvents A and B. To make contact with measurable quantities, 
it is convenient to express Vb(q) in terms of the Flory-Huggins interaction parameter XyEZl 

m - Va - (13) 

Combining equations (12) and (13) yields 

v (a) = S A 1 (q)S^ 1 (q) - 2 XA p5B 1 (q) ~ 2XBp5 A 1 (q) + D (u) 
S A - 1 (q)+Sv 1 (<l)-2xAB 

where D is a function of the corresponding Flory-Huggins interactionEI parameters only 

D = 2XABXAP + 2xabXbp + 2xapXbp - Xap - Xbp - Xab ( 15 ) 

Equation (14) could be written in a different form which is appealing from the point of view of the effective potential 
acting on a polymer chain in mixed solvents 

rrt \ 1 o ' 1 o (gg^g) ~ Xab - Xbp + Xap) 2 , . 

^°( q ) = ~Q~77t\ ~ 2 ^BP = -q-TT - 2 XBP - ! ^ " (16) 

^s(q) S B (q) S' A i (q) + 5 B x (q)-2xAB 

Similar results were obtained by others and in particular by Schultz and Flory§ and by Benoit and Strazielle@ 
using different methods of the mean field theory. The Edwards Hamiltonian method gives another route for deriving 
the effective potential on a chain together with the scattering properties. It has the merit of showing the response of a 
polymer chain undergoing interaction from the mixed solvent under various conditions of temperature, composition, 
etc. 



C. The scattering functions 

In the following subsection the crudest approximation that can be imagined. Although we consider only one chain 
in the mixed solution we use a mean field description of the RPA form. We are aware that this approximation is 
very bad indeed, nevertheless it is used here to sec the principal influence of the binary solvent on the interaction 
potentials. Has we used Gaussian chains this approximation would be fine. This case belongs to Vpp = 0. Another 
possibility of slight improvement is to treat the polymer plus excluded volume interaction as effectively " bare" system 
and use an appropriate Pade approximation including the selfavoiding walk exponent v « 3/5. 

Starting from equation 11, after some straightforward manipulations, one could write the Edwards Hamiltonian 
in terms of the scattering functions for incompressible mixtures of polymers and mixed solvents. The following 
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conventional notation is adopted. For interaction potentials, we keep the subscript P for all the polymer chains but 
for the single chain structure factor we use the subscript 1 meaning a single test chain. With this convention, the 
Edwards Hamiltonian becomes 



m = Y. K 1 ^) + V ™ + y BB + ^(q) - 2V BP ] \p\\ 



q 



+ [V A p - V BP - Vkb + ^bb + 5 B 1 (q)] p\p K 



qF-q 



+\ [Vaa + S^(d) + Vbb + 5 B 1 (q) - 2V AB ] \p$\ 2 (17) 
The above equation takes a greatly simplified form in the matrix notation 

m = \Y.p T A*r 1 p-* ( 18 ) 

q 

where p is a column vector 



is its transpose and 5(q) is a second rank square matrix whose elements depend upon the single particle structure 
factors 5 A (q), 5 B (q), 5°(q) and the Flory-Huggins interaction parameters Xab, Xbp, an d Xap- 

c-i (n \ = { 5f 1 (q) + 5 B 1 (q)-2x B p 5 B *(q) - xab - Xbp + Xap \ 

[q) Ub^^-Xab-xbp+xap 5 A 1 (q) + 5 B 1 (q)-2xAB J { ' 

Inversion of this matrix gives the partial structure factors 5n(q), Saa(c[) and 5Ai(q) 

g (q)= 5 A 1 (q) + 5 B 1 (q)-2xAB m) 

(Sf^q) + 5 B " 1 (q) - 2 X BP)(5 A 1 (q) + 5 B " 1 (q) - 2 X ab) - (5 B 1 (q) - xbp - xab + Xap) 2 

c (a) = 5r 1 (q) + 5 B 1 (q)-2xBP , . 

AM ' (5r 1 (q) + 5 B 1 (q)-2xBp)(5 A 1 (q) + 5 B 1 (q)-2xAB)-(5 B 1 (q)-XBP-XAB+XAp) 2 1 ' 



5 A i(q) = 5 1A (q) 

5b *(q) - xab - XBP + XAP 



(S^(q) + S B \q) - 2XBP)(5 A - I (q) + 5 B 1 (q) - 2 X Ab) - (S^q) - XBP - XAB + Xap) 2 



(23) 



These results can be used to extract the properties of the mixture both in the thermodynamic limit by letting 
q = and at finite q where we can explore the chain conformation and the spatial distribution of particles within 
the medium. We first consider the thermodynamic limit in the following section. 



1. The thermodynamic limit q = 

In this limit q — and equations (21) and (22) become 

e-i ,-i . ,-i ( , (0b 1 - Xab - Xbp + Xap) 2 \ 

Sl1 =01 + ^ - 2 { XBP+ 2(^ + ^-2xab) ) (24) 

q -l _ ,-1 ,_i , (^B 1 -XAB -XBP +XAp) 2 'j , 

+ ^ (^Xab + 2( ^ 1+ ^_ 2xBp) ) (25) 
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For simplicity, we remove the argument q from the structure factor in the thermodynamic limit q = 0. We consider 
a mixed solvent in the vicinity of the critical temperature T c where xab is close to the critical value Xc an d show 
how a small amount of polymer added to the mixture shifts the criticality condition. We also examine the way in 
which the chain responds to the approach towards solvent criticality condition. 
In the pure solvent mixture when q = 0, one has 

S AA = 2{ Xc - xab) (26) 

which can be obtained from equation (25) by letting cf>i = 0; Xc is the critical parameter for phase separation which 
is function of the solvent composition only 

Xc = I^a'+^b 1 ) (27) 

If polymer is added to the solution, the partial structure factor in the forward direction Saa becomes function of the 
polymer volume fraction <f>\ and can be written as 

S A l = 2( X ' C - xab) (28) 

This result shows that the critical parameter undergoes a shift by an amount Ax c 



Xab — Xbp + Xap) 2 

2(0r i A r r 1 + 0b 1 - 2xbp) 



Axc = X c - Xc = — ,"_i Ar _i ; ,_i — — r ( 29 ) 



One observes that the critical parameter x c = Xc + Ax c either increases or decreases depending upon the sign of 
Axc which in turn is fixed by the sign of the denominator in equation (29). This sign is determined by the relative 
magnitude of xbp as compared to 2 4>Ini 25^' ^ Xbp ^ ne interaction between polymer and solvent B is smaller 
than 24>In! ^ ne quantity Ax c is negative and the critical parameter decreases (i.e.x' c < Xc)- This means 

that the polymer favors solvent demixing since a smaller interaction drives the solvents to the limit of stability. 
For example, choosing N±= 100, <j>\= 0.01 and 0b=O.495, gives 2 4>\n 1 2^" = This means that if xbp is 

smaller than 1.5, Ax c is negative and addition of a small amount of polymer (approximately 1 percent) induces 
phase separation of the solvent mixture. If xbp is higher than 1.5, the system behaves differently since addition 
of polymer would result into a positive Ax c indicating compatibility enhancement of the solvent mixture. The 
mixture phase separates only under a strong repulsion between A and B solvents. These features can be observed in 
figure 1 where x'c is plotted against 4>\ for two values of xbp- One observes different tendencies, depending upon 
the value of the interaction parameter xbp- For the lower value, there is a slight decrease of x' c upon addition of 
polymer. If xbp is as high as 1.5 expressing a large polymer-solvent B incompatibility, a small amount of polymer 
would result into a large increase of Ax c meaning that the solvents A and B prefer to remain mixed rather than 
adsorbing on the polymer even though the interaction xab chosen in plotting this figure (xab = 1) is relatively 
high. This compatibility enhancement of the mixed solvent due to the presence of polymer is quite spectacular since 
Xc increases by practically an order of magnitude upon addition of a small amount of polymer. Qualitatively, this 
behavior is consistent with the observation made on PS/(benzene+cyclohexane) by Varra and AutalikO. These 
authors found that a strong interactions between PS and non solvents (cyclohexane) could lead to a reduced solvent- 
solvent interaction and hence a compatibility enhancement of the mixture. Furthermore, one finds that the shift 
x' c — Xc — Axc is only moderately sensitive to the interaction asymmetry between the polymer and the solvents 
A and B which is expressed by the parameter e = xbp — Xap- The shift in the critical parameter Ax c is strong 
upon addition of polymer even for small values of e. The interaction asymmetry parameter plays a more crucial role 
in connection with the polymer conformation near the solvent criticality condition as we shall see in the following 
section. 

The effective potential Vo of equation (16) is directly related to the classical excluded volume parameter acting 
on a polymer in the infinite dilute limit. This is made clear by writing the partial structure factor 5n(q) in the 
reciprocal form similar to Zimm's equation 

ST^q) = Si _1 (q) + (30) 
Letting Axo = Xc ~ Xab one obtains 
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VO - T ^XBP r-r (31) 

0b 2Axo 

For an upper critical solution temperature mixture, the interaction parameter xab increases approaching \c- In this 
case, the polymer reacts in different ways depending essentially on whether e is zero or not. This is a clear illustra- 
tion of the importance of interaction asymmetry between polymer and solvents in determining the thermodynamic 
behavior of the solution. As the solvent mixture approaches the criticality condition (Axo — * 0), a striking difference 
exists between the case where the solvents present the same interaction with the polymer (e = 0) and the case where 
the two solvents interact differently with the polymer (e ^ 0). Equation (31) shows that in the case of a symmetric 
interaction (e = 0), Vo decreases linearly with Axo- However, in the non symmetric interaction problem (e ^ 0), the 
third term in the RHS of equation (31) increases rapidly as Axo — > and Vo undergoes a rapid decrease. Following 
this sharp decrease, the effective potential vanishes and may even become negative as xab approaches Xc- The chain 
undergoes a phase transition from a swollen state to unperturbed dimensions or even a collapsed configuration. This 
is illustrated in figure 2 where Vo is represented as a function of Axo for e = and 0.2. The upper curve in this 
figure corresponds to the case where the interaction between polymer and solvents A and B are the same and e = 0. 
It shows that Vo increases relatively smoothly when Axo — » 0. 

One could also introduce an effective polymer-solvent interaction parameter xo via 

v o = -t 2 Xo (32) 

<PB 

With 

(Axo + e) 2 . . 

Xo = Xbp + 4Axp (33) 

Figure 3 represents the variation of the effective polymer-solvent interaction parameter xo as a function of Axo for 
the same values of e as in figure 2 and shows similar tendencies. For e = 0, xo decreases slightly when Axo - * 0, 
while for e = 0.2, there is an inversion in the variation of xo expressed by a turn up and a sharp increase as the 
mixed solvent approaches criticality indicating a strong polymer-solvent repulsion. 



2. The small q limit 

As a first approximation, we consider that q^ 1 ^> 1 which means that by scattering radiation, one observes the 
solvent particles as point like and the q- dependence comes entirely from the form factor of the polymer only. This 
approximation could be easily relaxed if one wants to include intramolecular interferences within the solvent particles 
themselves and explore by radiation scattering the internal structure of chains. Moreover, we will assume that the 
form factor of the chain could be approximated by 

P(q) = — W (34) 

1 + ^ 

Following these assumptions, equation (17) can be written in the form 
Sii(g = ) 

i + q 2 e 



= (35) 



Where the correlation length £o represents the distance over which the polymer configuration is sensitive to solvent 
fluctuations. It depends upon the polymer radius of gyration R g and the excluded volume interaction Vo according 
to the relationship 

R 2 

e ° = 2(1 + Wi) (36) 

Figure(4) represents the variation of £o/Rg as a function of Axo for e = and 0.2. For e = 0, S,o/Rg decreases 
continuously when Axo — » while for e = 0.2, £,o/Rg increases suddenly when Axo approaches zero. One can 
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visualize the polymer as a chain of blobs of size £o- Over distances below the blob size, the chain is sensitive to 
solvent fluctuations and its size reduces as a result of increasing effective interaction xo- Over distances exceeding the 
blob size, the chain does not feel the effects of solvent fluctuations and shows a comparatively swollen conformation. 
As the solvent approaches criticality when A\o — > 0, the correlation length £o becomes comparable or larger than 
R g implying that the entire chain is subject to solvent fluctuations and undergoes a phase transition from a swollen 
to a more compact conformation. For e = and xbp = Xap = 0.5, the chain feels theta solvent conditions even 
though Axo — > and solvents demixing conditions are approached. n 

The blob picture for a polymer in mixed solvents was first suggested by Brochard and de GennesH who examined 
the polymer conformation in a mixture of two good solvents, when the polymer affinity differs substantially for the 
two solvents. They were the first to predict that the polymer adopts a collapsed configuration near T c even it would 
be swollen in either of the pure solvents under similar conditions. 

3. Many chains in mixed solvents 

A. The effective potential of a test chain 

In this section, we discuss the problem of many chains by extending the Edwards Hamiltonian formalism to strong 
solutions of polymers and mixed solvents. The effective potential acting on a test chain together with partial 
structure factors are calculated as a function of polymer concentration and mixed solvent properties. The Hamiltonian 
formalism is extended by taking into account the effects of the many chains present in the medium. This gives the 
following result for the partition function including summations over all chains designated by the letters a and (3 

"o^E / / 5{K a {s)-TLp{s'))dsds' 

-^apE/ E^w-^-^E / E 5 ( R «( s )- r ?) ds 
-\vaa E^ - 4) \vbb E *(>f - 4) - ^ B E - r f )] ( 37 ) 

ij ij ij 

The subscript a and (3 run over all chains in the solution except a test chain designated by the subscript 1 (meaning 
single). Transforming space coordinates into collective density variables leads to the partition function 

-V P P E IPqllPql - \ E^PP + ^(COKI 2 

q q 

-v AF E p>\ vap E pIp\ - y** E p>\ - E \pI\\pI\ 

q q q q 

A E^AA + ^ X (q))|pt| 2 - \ E^BB + 5 , B 1 (q))|Pq| 2 - E *W£p? q ( 38 ) 



2 ^ ™ v "" lrql 2 

q q 



Where p A anc i pB^ ^p anc j ^ are ^ ne collectives coordinates for the solvents A and B, the polymer P and the test 
chain f, respectively. In order to identify the effective potential acting on a test chain together with the partial 
scattering functions for the multicomponent mixture under investigation, it is convenient to write the partition 
function in matrix form 
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£=/m l(S )/n^exp[-j^ (^) ^->pekh 

exp[- £ 2V(q)p q - J2 ^M(q)p- q ] (39) 



where /? q and V(q) represent the column vectors 




(40) 



jVap J* exp(-iqR(s)ds) 



2 



V(q)= | iV B pJ^e X p(-iqR(s)ds) | (41) 
IVpp / exp(-iqR(s)ds) 

and M(q) is a square three-by-three matrix 

/ K^aa + ^j) |Vab 
M(q)= ^ab UVbb+s^) ^bp j (42) 

V 5^AP ^ 

Integrating over the collective variables p q , one obtains the partition function as follows 

Z = J PR l(s) ex P [-JL^ (^j^-^-E^^ + E^M^ 
Recalling that Z = J T)Ri(s) exp— f3H yields 

£ (^7^) 2 ds + E ^ ^ ^ ^P (~iq(R(s) - R(s'))) (44) 



(43) 



2/ 2 

where one can identify the effective potential Vq acting on a test chain. If the volumes occupied by monomers and 
solvent molecules are the same, then Vpp, Vaa and Vbb are approximately represented by the same quantity denoted 
Voo and one can introduce the Flory-Huggins interaction parameters using Vqo — Vij - Xij where the subscripts i and 
j run over A, B and P. All polymer chains in the medium have the same interaction parameters and no distinction 
exists between the test chain 1 and the others. The effective potential V(q) is not only due to the solvents but also to 
all the chains present in the medium. For incompressible mixtures, the potential Voo goes to infinity and one obtains 



v , , = ^(q^B^q) + d- 2xApSB 1 (q) - gxgp^q) (45) 

' 11 5A 1 (q)+5B 1 (q)-2xAB + *(q)(5 A 1 (q)^ B 1 (q)+D-2xAp5 B 1 (q)-2XBP^ A 1 (q)) 



Here D is defined in equation (15). It is interesting to note that the effective potentials acting on chains in infinitely 
dilute solutions Vb(q) and in strong solutions V(q) are related via 

V-^q) =y - 1 (q)+5 P (q) (46) 

This result is characteristic of the series expansion in the Ornstein-Zernike modelEl It suggests that one can write 
the total potential V(q) as an infinite series of interaction terms. The order of these terms depends on the number 
of intermediate chains involved in the interaction (see figure 5a). Each diagram corresponds to a particular term in 
the series expansion 

V(q) = ^o(q) - Vo(q)S P (q)Vo(q) + yo(q)Sp(q)^o(q)5 P (q)V r o(q) + ... (47) 
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These series can be summed up in a similar way as in the Ornstein-Zernike expansion ( see figure 5b) . The summation 
gives 

V(q) = Vb(q) - V (q)S P (q)V(q) (48) 
which is equivalent to 

V(q) = V (q) [1 + F (q)^p(q)] _1 (49) 
B. The Edwards screening 

For completeness, we remark that the Edwards screening can be recovered. The effective potential for many chains 
in equation (47) can also be written in the following form 



V(q) = V (q) - ;° KH [ r/ - (50) 



^ 2 (q) 

S p l (q) + V Q (q) 

Using Edwards prescription for the form factor of Gaussian chains 
P JVp 

Sp = TTWiq (51) 

where the radius of gyration is related to the degree of polymerization N and segment length / by the known result 
i? g = ly/N/6. Substituting equation 51 into 50 yields 

V q 2 t? 

V Q = TTg^a (52) 
Where £ is the Edwards screening length 

v / 12Vo(/)p 

In the infinite dilute limit, <fip goes to zero and the screening length £ tends to infinity. As one would expect, the 
effective potential V tends to Vq in this limit. Figure 6 represents the ratio V/Vq versus q£ as given by equation 52. 
The ratio V/Vo is small for g£ less than 1 indicating screening of the potential field acting on a test chain. When q£ 
increases the ratio V/Vq tends to 1 and the effective potential is practically the same as in infinite dilute solutions. 
Locally at large q, the test chain behaves as if it were isolated and is essentially not perturbed by the presence of 
other chains in the medium. It feels only the fluctuations of the solvent mixture in its immediate vicinity. 

This analysis of screening presented here and following the Edwards prescription is only qualitative. A more precise 
treatment is needed if one wants to have a quantitative evaluation of the screening effects. Such an impressment can 
be made using other methods such as renormalization group theory methods or field theoretical tool: 



BBS" 



C. The scattering functions 

The scattering matrix can be readily extracted from the partition function in the matrix notation 

/ 5r 1 (q) + 5 B 1 (q)-2xBP 5 B 1 (q)-2 X BP S B 1 (q) - X AB - XBP + XAP \ 

S" 1 (q)= 5 B 1 (q)-2 XB p 5 p 1 (q) + 5 B 1 (q)-2 XB p (q) - Xab - Xbp + Xap (54) 

V ^B 1 ^) -XAB -XBP +XAP S B (q) - XAB - XBP + XAP 1 (0> + ^B 1 (<l) ~ 2 XAB / 

The diagonal elements of the scattering matrix S(q) define the contribution to the scattering signal due to dif- 
ferent constituents in the mixture. Three quantities are particularly relevant for our purpose in the present work. 
Considering the solvent B as a background medium, one can obtain the diagonal elements of equation (54) denoted 
Sii(q), Spp(q) and S < AA(q). The properties of these partial structure factors are briefly discussed below. 



10 



1. Structure factor of a test chain Sii(q) 



The first element of the matrix in equation (54) gives the scattering function from the test chain 

S 1 -i 1 (q)=Sr 1 (q)+^i(q) (55) 

where Vi(q) is identical to V(q) given in equation (46). The subscript 1 is introduced to describe a generalized 
excluded volume experienced by a test chain labeled 1 in the strong solution. Discussion of the screening presented 
above could be applied here invoking the excluded volume interaction experienced by the test chain 1. It is worthwhile 
to note that if the mixed solvent satisfies the criticality condition in the absence of polymer, then + = 2xab, 
and equations (55) and (45) give 

1 1 1 

~- = it + c- (56) 
ill bp 

This result corresponds to the scattering function, one obtains from a solution of identical chains under theta 
conditions where excluded volume interactions are completely screened out. It is consistent with the variation of Vq 
versus A% which shows that the apparent excluded volume tends to as A% — * {figure 2) 



2. Structure factor of the polymer Spp(q) 

The second diagonal term of the scattering matrix 5(q) represents the contribution from the polymer and found as 

Sp^(q) = S p 1 (q) + V P (q) (57) 

Where Vp(q) has a similar expression as Vi(q) except that Sp(q) and Si(q) are interchanged. If the mixed solvent 
approaches criticality, the signal from the p-chains become identical to the signal from the test chain 1 and hence 
the same as in equation (56). The chains are under theta condition as indicated earlier. 



3. Structure factor of the mixed solvent Saa(q) 

The third diagonal element of the matrix in equation (54) gives the partial structure factor of the mixed solvent 

AAW BW XAB (^p 1 (q) + 5r 1 (q))(5B 1 (q)-2XBp) + 5r 1 (q)5p 1 (q) k ' 

In the infinite dilute limit Sp = and one recovers the former result given in equation (22). Considering the q 
= limit, one can make similar remarks as those following equation (29), except that the concentration range of 
polymer is not limited to infinite dilute solutions and the polymer concentration range is extended to higher values 
in investigating the variation of the interaction parameter x' c as a function of polymer concentration <f>p. 



4. Conclusions 



In this paper, we present some theoretical aspects of polymers in mixed solvents using the Edwards Hamiltonian 
formalism. Discussions are confined to thermodynamic and structural properties and both the single and the many 
chain problems are considered. If one solvent species presents a strong repulsion towards the polymer, addition of 
a small amount of polymer leads to compatibility enhancement of the solvent mixture. The critical parameter at 
which the solvent mixture phase separates increases dramatically. The effective potential acting on the polymer 
in the medium depends enormously on the asymmetry in the affinity of the polymer with respect to both solvent 
species. If this asymmetry is large, as the mixed solvent approaches the criticality condition, the effective interaction 
polymer-solvent increases dramatically. The polymer undergoes a phase transition from a swollen configuration 
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to unperturbed dimensions and eventually to a collapsed configuration as the mixed solvent approaches criticality. 
When the correlation length for solvent fluctuations is of the order or higher than the radius of gyration of polymer 
R g , the chain conformation becomes characteristic of a polymer in bad solvents and the effective excluded volume 
parameter vanishes and may even become negative. 

In the many chain problem of strong solutions, the effective potential acting on a test chain is calculated using the 
same formalism based on the Edwards Hamiltonian. This potential is related to the bare effective potential in the 
infinite dilute limit via a series expansion which is reminiscent of the Ornstcin-Zernike model. This series expansion 
is represented graphically via a diagrammatic method where the successive diagrams represent interactions involving 
increasing numbers of intermediate chains. If the solvent mixture is critical, interesting results are obtained for 
the structure factor of the polymer and a test chain. If the solvent mixture is near criticality, all excluded volume 
interactions are screened out and chains in the medium feel theta solvent conditions. 

This investigation is based upon the Edwards Hamiltonian formalism and gives a qualitative picture of the thermo- 
dynamic behavior and structural properties of polymers in mixed solvents. Near critical regions where the fluctuations 
are strong, one has to rely on an appropriate model to achieve a more rigorous quantitative evaluation of the effects 
of critical fluctuations on the thermodynamic and the scattering properties of polymers in mixed solvents. This could 
be done within the framework of field theoretical or rcnormalization group theory calculations. 
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Figure Captions 

Figure 1: 

Variation of x' c the critical interaction parameter between solvents A and B with the polymer volume fraction <f>i for 
two values of xbp and xab = 1, iV = 100, 4>b — 0.495 e = xap — Xbp = 0.2 (see equation (29)). 

Figure 2: 

Variation of V the effective excluded volume parameter with the parameter Axo for two values of e and Xbp =0.5, 
(see equation (31), the other values are chosen to be the same as above) 

Figure 3: 

Variation of xo the effective interaction parameter with the parameter Axo for two values of e and xbp = 0.5, (see 
equation 33) 

Figure 4: Variation of £o/-Rg the normalized correlation length for solvent fluctuations with the parameter A%o for 
two values of e and xbp = 0.5, (see equation (36)) In plotting these figures, we have used N = 100, 0a = 4>b and 

XAB = I- 

Figure 5: Diagrammatic representation of the relationship between the effective potentials in the many chain 
(V) and in the single chain problems (Vb) a) Representation of the infinite series of equation (47) b) Representation 
of equation (48). Vertical thick lines represent the main interacting chain. Vertical dashed lines represent intermedi- 
ate chains of interaction. Horizontal dashed wavy line is an interaction Vb Horizontal thick wavy line is an interaction 
V. 

Figure 6: Variation of the ratio V/Vq as a function of q£ according to equation (52) where £ is the Edwards 
screening length. 
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